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Abstract 



We consider various ways of treating the infrared divergence which appears in the 
dynamically generated fermion mass, when the transverse part of the photon propagator 
in N flavour QED3 at finite temperature is included in the Matsubara formalism. This 
divergence is likely to be an artefact of taking into account only the leading order term 
in the expansion when we calculate the photon propagator and is handled here phe- 



nomenologically by means of an infrared cutoff. Inserting both the longitudinal and the 
transverse part of the photon propagator in the Schwinger-Dyson equation we find the 
dependence of the dynamically generated fermion mass on the temperature and the cutoff 
parameters. It turns out that consistency with certain statistical physics arguments im- 
poses conditions on the cutoff parameters. For parameters in the allowed range of values 
we find that the ratio r = 2 * Mass(T = 0)/ critical temperature is approximately 6, 
consistent with previous calculations which neglected the transverse photon contribution. 
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1 Introduction 



At zero temperature, a considerable quantity of work has been done |l -HJ and has been 
extended to finite temperature more recently |14|-|T9|j. Our interest here is to carry forward 
the finite temperature calculations by improving upon the approximations made, as described 
below. In particular, we shall attempt to include the transverse photon modes, which have so 
far been neglected in all the calculations done using the Matsubara formalism. 

At finite temperature Euclidean QED 3 in the reducible representation, there are three sig- 
nificant parameters T c , N c and r which need to be calculated. T c and N c are respectively the 
temperature and the number of flavours above which chiral symmetry is restored and at which 
there is a phase transition from a superconducting phase to a normal phase, r is defined as 
the ratio of twice the zero temperature fermion mass to the critical temperature, T c . This 
ratio is characteristic of the mechanism which leads to superconductivity. In Dorey and Mavro- 



matos [f[4|L5l a calculation was done using the Schwinger-Dyson equations at finite temperature 
and the Matsubara formalism. In the following we shall be using the term photon to describe 
not the electromagnetic photon but the statistical photon. The approximations made in \TA ) I3| 
are as follows: 

(i) The full photon propagator A MV was calculated to leading order in in Landau gauge. 

A m (|p|,/?) = ^ + ^ (1) 

where A^ u and B^„, the longitudinal and transverse projection operators, are given by 

a f x P»Po\ (p 2 \ ( . PoP A 

B, v = 6* - 8 ju (2) 

and IT™, II m are related to the one loop photon self-energy as follows: II™ = IT^A^ + IT^B^. 

(ii) Only the \i — v — part of A™ u was retained. 

(iii) Only the m = (pq = 0) part of was retained, retardation effects being neglected. 

Approximations (ii) and (iii) will be referred to as the instantaneous exchange approxima- 
tion which corresponds to considering the quantum statistical mechanics of an ensemble of 
Dirac fermions interacting via the static two-body potential instead of a full field-theoretical 
treatment. 
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(iv) The fermion mass E was assumed to depend only on temperature and not on frequency 
or momentum. 

£ w (|p|,/3)«£°(0,/3) = E(/3). 

(v) the fermion wave function A was neglected, A — 1. 

As a consequence of (iv), (v) and the Ward-Takahashi identity, the full vertex T M was taken 
to be e7 M . 

The results were that there is dynamically generated mass which is a decreasing function 
of the temperature T and vanishes for temperatures above a critical temperature T > T c , T c 
is itself a decreasing function of the number of flavours N. No critical number of flavours 
N c was found when the zero temperature limit was taken and no N c was expected at finite 
temperature. This is in agreement with the zero temperature calculation of [|IJ, because this 
work adopts the constant mass approximation (iv) as well. The ratio r was found to fluctuate 
around 10 as N increases. This is a large value compared to the BCS r ~ 3.54 and seems to 
be phenomenologically correct, as there are indications that the layered CuO superconductors 
indeed have larger r than BCS systems. The origin of this phenomenon is traced to the 
appearance of a temperature dependent mass for the temporal photon (plasmon mass), which 
implies that external electric fields are screened by thermal excitations. 

In Aitchison, Dorey, Klein-Kreisler and Mavromatos |L6], the constant mass approximation 



(iv) was relaxed and the momentum dependence was taken into account 

E m (|p|,/?)«E (|p|,/3). 



The results were similar to those of [p~4| , |15|] , the only difference being that a temperature 
dependent critical number of flavours N C (T) was found. This is the same behaviour seen 
at zero temperature when momentum dependence was taken into account [|. The zero 
temperature limit of the N C (T) found in |1| is just greater than 2 and therefore does not 
coincide with the value 3.2 of [@||. This is due to the fact that the finite temperature calculation 
of (HJ employed the instantaneous exchange approximation which was not used in the zero 
temperature calculations p,|3|] . It should be stressed that the relaxation of (iv) did not affect 
the value of r which remained around 10. 

Following the zero temperature arguments of 0-0], Aitchison and Klein-Kreisler [|I7| per- 



formed a finite temperature calculation dropping approximations (iv) and (v) simultaneously. 
The Pennington- Webb vertex T v = ^ V A of j|,[5| was adopted. Unlike the more elaborate vertex 
of [|7] (Ball-Chiu vertex + a non-trivial transverse part), the Pennington- Webb vertex does 
not satisfy the Ward-Takahashi and Ward identities. However it yields results which are in 
qualitative agreement with those of 0. Although works show that taking into account 
the momentum dependence of the mass and the wave-function renormalisation leads to the 



conclusion that there is no N c at zero temperature, |17| suggests that this not the case at fi- 
nite temperature. The reason can be traced to the different infrared behaviour of the photon 
self-energy at zero and at finite temperature. Thus the results of [17] are very similar to those 
of the earlier finite temperature work |TB| and the value of the ratio r was once again found to 
be approximately 10. 

The r = 10 feature of QED3 seemed to be remarkably robust, surviving all refinements, 



until Aitchison |18[ showed that dropping the instantaneous exchange approximation (ii), (iii) 
leads to a considerable reduction of the value of r from 10 to 6. One obvious problem with 
the instantaneous approximation is that it cannot reproduce the well-studied zero tempera- 
ture limit. This is because as T — > 0, all frequency components are to be included; but the 
instantaneous approximation retains only the n = component of the photon propagator, so 
the solution cannot join continuously on to the T = solution. A second problem (indicated 
here for the first time) may be that the instantaneous approximation cannot reproduce the 
requirement that ^ — > as T — ► 0; we shall show in Section 5 that this requirement must be 
true of physical solutions, if the thermodynamic arguments of that section hold. 

The first study of retardation effects in the S-D equation in QED3 at finite temperature 



18| , was performed in the real time formalism. This has the advantage that there is no 



discrete sum over frequency components, so that the zero temperature result could be recovered 
straightforwardly. The reduced value 6 was calculated for r. No critical number of flavours N c 
was found. This calculation, however, used an approximate form for the photon propagator, 
which did not treat correctly the non-analytic structure near the origin. Therefore some doubt 
remains as regards the reliability of the results. 

Finally Lee in made an attempt to drop the instantaneous exchange approximation (ii) 
and (iii) in the Matsubara formalism this time. Convenient approximate forms for both the 
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longitudinal and the transverse part of the full photon propagator were found to leading order 
in 1/N, having the correct behaviour for small 3-momentum. Both the constant mass approx- 
imation and a second approximation, where the mass was taken to be frequency dependent, 
were investigated. As expected, in the former case no N c was found whereas in the latter case 



there was a N c . In both cases r ~ 6 in agreement with |18fl . However the problem with |L9[ is 
that, in the S-D equation for the full fermion propagator, the transverse part of the full photon 
propagator was discarded, because of its divergent infrared behaviour. This new approxima- 
tion is qualitatively similar to (ii), since has no transverse part as well. Thus this paper, 
although it claims to discard both (ii) and (iii), in reality drops only (iii) and substitutes (ii) 
with an equivalent approximation. 

Our purpose is to go beyond |1| by including the transverse part of the full photon propa- 
gator in the S-D equation for the full fermion propagator at finite temperature. In other words 
we shall completely discard (ii) and (iii) but retain (i), (iv) and (v). 

The infrared divergence of the transverse contribution of the full photon propagator makes 
the S-D equation meaningless and we must find some way of dealing with it. One might at 
first think that the trouble has arisen, because the photon self-energy has been calculated with 
a massless fermion in the loop, instead of the self-consistently generated mass E. However 
including E in the loop calculation does not help, as we show in Appendix A. We are therefore 
forced to consider the introduction of some form of infrared cutoff, if we want to proceed at 
all. We shall consider two types of cutoff: the first is on the momentum integral of the S-D 
equation, which we place at |k| = e/c#T; the second is a mass cutoff in which we give the zeroth 
mode of the transverse photon contribution a squared mass of the form Saln2kBT similar to 
the plasmon mass of the temporal part of the photon propagator. We may view these cutoffs 
as a phenomenological way of taking into account terms beyond the leading order in 1/N, since 
it seems likely that the infrared divergence would be removed, if the full (all orders in 1/N) 
zero frequency contribution could be calculated. The introduction of infrared cutoffs of both 
the integral and the mass forms is a usual practice in the zero temperature literature, although 
there is no infrared divergence in this case Our purpose then, is to investigate how the 

fermion mass E depends on the two types of cutoff and to attempt to extract sensible physical 
conclusions from these calculations. 
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In Section 2 we formulate the gap equation for E. In Section 3 we consider what we call 
"extreme" solutions in which the infrared cutoffs e and 5 are taken to be so small that the 
infrared divergence is the dominant effect. In Section 4 we present numerical solutions for a 
range of values of e and 5. We find that there are two classes of solution: those at "large" e 
and 5 for which < for all T and those at "small" e and 5 for which ^ > at small T. 
In Section 5 we shall argue that the latter class is unphysical by considering the entropy of a 
dilute gas of fermions with mass S(T). Accepting this argument, it follows that there exists 
critical values e c and 5 C below which the solutions are unphysical, in the sense of ^ > for 
small T. These critical values are obtained in Section 6 by finding a simple analytical form 
for the numerically calculated E(T, e or 5), which may itself be useful in other contexts. Some 
final comments are made in Section 7. 



2 The Schwinger-Dyson equation for E 

The Lagrangian of massless, Euclidean QED 3 with N fermion flavours is 
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L = --f^r + J2$i(i0-em (3) 
4 1=1 

We consider the reducible representation of the Dirac algebra so that the above Lagrangian 
has two continuous chiral symmetries. In the Matsubara formalism, the full Schwinger-Dyson 
equation for the full fermion propagator at finite temperature takes the form 

s F \p f ) = sf-\ Pf ) -- E / ^i„s F (k f )AMr»(k f , Pf ). (4) 

where we have used the subscripts / and b to denote fermionic and bosonic momenta respectively 
and 



(2n + l)7r 

Pf = (Po, P) Po = p 

(2m + 1)tt 



k f = (k , k) k 
q b =p f -k f = (g , q) go 





2{n — m)-n 
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After applying approximations (i), (iv), (v), taking the spinor space trace and performing the 
trivial angular integration, (^) reads 



where 



S(/3) 



N(3 ^ 



£(/?) 



2tt 



(2to+1)tt 

P 



+ k 2 + E(^) 1 



(5) 



MM 



(|k|,/3) 



1 



2m7r 



+ 



1 



2m,7r 



k 2 + n^(|k|,/?) 



(6) 



+ k2 + nr(|k|,/3) 

A trivial solution of (j5|) is that £ is identically zero for any temperature. However this solution 
for the fermion mass would rule out QED3 as a model of superconductivity. If QED3 is to be 
such a model, it should be possible to prove that there exists a critical temperature T c (see 
Appendix B) below which there is a non-zero solution of (|5|). Thus we proceed by assuming 
that £ is not identically zero and seeking a self-consistent solution of the resulting equation: 



a 



E 



|k|d|k| 



A-"V(|k|,/% 



1 



(2to+1)tt 

P 



k 2 + E(P)'< 



(7) 



Both (|7|) and are well-known and appear for example in [14j, 15]. The exact expressions 
for n^ m and n^ m are also given in \T%\. However (|?p is too hard to attempt even a numerical 
solution. This is exactly why a further approximation, such as the instantaneous exchange 
approximation, was necessary. Recently in |19[ it was shown that a less crude approximation 
is sufficient to render (|7|) numerically solvable. 



For m = 



n° 



n? 



a \k\(3 



. lo 16 In 2 

K\/3 H e 16 la 2 



|k|/3 



7T 



g 16 In 2 



For m/0 



a 



A 



2m7r 



/3 



k 2 . 



As noted in p|, when /3|p| -> 0, n° -> /^ 2 |pp|. / is a constant which is found to be n 
in our approximation of the photon propagator; this differs by a factor of 3/4 from the |p| — » 



behaviour of the exact result (see Appendix A). For n° the low |p| behaviour is different; in 
" hl2 §, this constant value being the plasmon mass. 



this limit, n° 



7T ft' 
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As we have already explained, [TS[] does not make full use of these more accurate forms and 
in particular of the expressions for n'g and IT§,m 7^ 0. The approximate form for 11^ had 



also been employed in the earlier work JT6[ in combination with the instantaneous exchange 

approximation. 

On rearrangement (H) reads 



1 



2A% 



[S L (a,s) + S T (a, s)]. 



(9) 



St{o>, s) is the contribution from the transverse part and takes the form: 



/■oo 

St(cl, s) — / x dx 
Jo 



1 1 



x 2 + /3 2 U° 1 x 2 + 7i 2 + a 2 s 2 
1 



X 



^ x 2 + (27rm) 2 + 0.125a[x 2 + {27rm) 2 ] 1 / 2 

v uo, 



y x 2 + [2n(m + 1/2)] 2 + a 2 s 2 x 2 + [2n(m - 1/2)] 2 + a 2 s 2 1 

where x = |k|/3, a = a(3 and s = E/a. 5^ is the longitudinal contribution to (0) and is exactly 
the same as St in ([10|), except for II? being replaced by II3. 

The low |k| behaviour of II^ is oc k 2 . It follows that there is an infrared divergence in 
([T0|), which is absent in the case of Sl due to the presence of the plasmon mass in II3. We 
note that using the finite fermion mass self-energy, 11^ (Appendix A), instead of the zero mass 
quantity IT^, would not cure this divergence. We are left with one other alternative at leading 
order in 1/N: to impose an infrared cutoff. 

We shall look at two types of cutoff. In the first, which we call the integral cutoff, we 
replace the lower limit of the |k|-integral in (|7|) by e/c^T, e being an arbitrary constant. On 
rearrangement of the terms leading to (|9p and ( [TDD 5 thi s cutoff becomes the lower limit of the 
a;-integral, e. The idea behind choosing such a cutoff is that at T = there is no infrared 
divergence in (^]) and we shall want our T = solution to join smoothly on to our T ^ 
solutions, when we consider the contribution from all the modes in Sl and St- The form eksT 
is the simplest cutoff on the |k| integral which ensures this, for as T — > the cutoff on the 
I k| -integral tends smoothly to zero. 
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In the second type of cutoff we shall add an extra plasmon-like term 5aln2/j3 to U®, where 
again 5 is an arbitrary constant; we call this procedure the mass cutoff. Again at T = this 
will vanish leaving a smooth T — > limit. 

We are able to evaluate analytically all the x-integrals for the m/fl modes of (§). In our 
evaluation of the x-integrals in the sums, we use an approximation: we neglect e in the integrals 
of all the modes except the m = mode in St and Sl- This approximation is only valid for 
e < 2tt, which is not a problem, as e is an infrared cutoff. 

Using this further approximation gives us the following expression in the integral cutoff case: 



_ Q(a,s) _ a ^ 
~ 2ttN ~ 2^iV X 



2/(27rm, 0.125a, (aV + [2vr(m + 1/2)] 2 ) 1/2 ) 

n=l 

+ 2/(27rm, 0.125a, (aV + [2vr(m - 1/2)] 2 ) 1/2 ) 

m=l 

i: 



r°° X dx 



x 2 + 7T 2 + a 2 s 2 \x 2 + n°(x)/3 2 x 2 + Ii Q 1 {x)(3 2 



[ID 



where 



I(d, a, c) 



2(a 2 + c 2 - rf 2 ) 



In 



(d + a) 



+ 



; c 2 _ a>2) 1/2 ( a 2 + c 2 _ ^2) 



arctan 



(12) 

For the mass cutoff our equation is similar but with the zeroth mode contribution in ([TTf) being 
replaced by: 



x dx 



/o x 2 + 7T 2 + a 2 s 2 (x 2 + ng(x)/3 2 + x 2 + n?(x)/5 2 + a5 In 2 J ' 
Although we have argued that terms of higher order in 1/N would regulate the infrared di- 
vergence, we have no way of telling whether the mass or integral cutoff is the best way of 
approximating such terms (although we shall see later that of these two cutoffs, if QED 3 is to 
be a model of superconductivity, the integral cutoff is preferred) or what values of e or 5 give 
the best approximation without explicitly calculating these terms. We can however start by 
considering the limiting case of very small e or 5, where we can make some progress analytically. 
We shall find in Section 3 what appear to be clearly unphysical features of the solution, thus 
confirming that the divergence is physically serious and suggesting that e (or 5 ) cannot be 
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too small for consistent physical results. This will lead us to a numerical exploration of the 
solutions for various e(S) in Section 4. We shall find that for e(5) less than a certain critical 
value e c (5 C ) the quantity ^ > as T — > 0, which is physically undesirable, as we shall explain 
by a thermodynamic argument in Section 5. Requiring < for all T therefore imposes a 
constraint on the possible smallness of e(S) and hence (on our interpretation) of the higher order 
contribution in 1/N. We shall also consider how sensitive the ratio r = is to variations of 
e(S) in the allowed regions e > e c , 5 > S c . 

3 The extreme solutions for E 

In this section we shall be interested in the behaviour of S for very small 5 and e for which we 
can see exactly what behaviour the infrared divergence induces in E. 

To find the solution of (|TT| ) at extremely small e(S) we need only retain the zeroth transverse 
mode; at sufficiently small e(8) this dominates fllT]) (at large enough Tks/a) so that all other 
terms may be neglected. We are therefore left with an equation - in the case of an integral 
cutoff - of the form 



Now we know that for small x, Tl\(x)f3 2 goes as a fx 2 . Since we are interested only in the 
singular contribution, when e — > 0, we may impose an ultraviolet cutoff A. The purpose of A is 
merely to simplify the analysis: below A we shall, to a good approximation, be able to replace 
with a fx 2 ; above A, we shall be able to neglect contributions from the various momenta, due 
to the fact that, for extremely small e (and 5), the integral from A to oo is negligible compared 
to the infrared divergent integral from e to A. By replacing oo by A we are able to write flli]) 
as: 




(14) 



2irN 1 + faJe x(x 2 + ir 2 + a 2 s 2 ) ' 
The integral in (^) can be easily evaluated, yielding the following equation for s: 




(15) 



1 




1 1 




(16) 
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We note that, since A is the ultraviolet and e the infrared cutoff, 



A A 2 + vr 2 + a 2 s 2 , A , A 2 + vr 2 + a?s 2 

A > e & - > W— — 4»ln- > In 

e V e + 7T + o s e 



e 2 + 7r 2 + a 2 s 2 



which means that the second term in can be neglected, leaving us with 



its solution being 



2nN(l + af) \n 2 + a 2 s 2 J ^ \e 



;i7) 



S(T) = T\ 



7T Z 



2irN(T + af) 
where h = In - . 

Later we shall comment on this solution but first we turn our attention to the mass cutoff 
case. Here the dominant contribution gives the equation 



2ttN Jo {x 2 + afx 2 + aS In 2) (x 2 + vr 2 + a 2 s 2 ) ^ 
where again we have imposed the ultraviolet cutoff A. By neglecting all terms which are not 
divergent, when 5 — > 0, we get 



In 



a8 ln2 + A 2 (l + a/) 
a51n2 



(20) 



4irN[(n 2 + a 2 s 2 )(l + af)-a5\n2] 

Now, if 5 <C A 2 / 2 and 5 <C which will certainly be the case for extremely small 5, we 
obtain the solution: 



E(T) 



T, 



\ 4ttA^(T + a/) 



In 



A 2 (r + a/) 
ar<Jln2 



7T Z 



T 



a 



A T 
In ; +lnW— + f 
\\ 2nN(T + af) [ v^hi2 V a J 



(21) 



Now we compare the two logarithms appearing in the second line of fl2~I|) . Dynamical mass 
generation (E ^ 0) occurs for T < T c whereas A is an ultraviolet cutoff. Furthermore a > e 2 
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whereas 5 is an infrared cutoff, f is just a constant. Thus it is reasonable to neglect the second 
logarithm. 



£( ^^ M^ -'' (22) 
where h = In . 

We see that fl22 ) is very similar to (|P8|) , fig.[l]b and fig.|l]a respectively. We can therefore 
conclude that, for very small cutoff and provided it vanishes at T = 0, the choice of cutoff type 
has little effect on the solution. In both cases the infrared divergence is present in S(T) as 
h or h — > oo. 

We note that (|18|), fl22|) define critical temperatures T c = — af, T c = — af 



2n 3 N > ±c 2 -k 3 N 

respectively. However E vanishes at T = as well. This is due to our having retained only 
the zeroth transverse mode: a finite value of £ at T = would be obtained, if we included the 
modes which are not infrared divergent. The real pathology of (pl|) or (|22|) (apart from the 
infinity in the limit as e or 5 — > 0) is that these solutions for S increase as T moves away from 
zero. This seems intuitively unnatural; the mass is analogous to an order parameter and we do 
not expect this to increase in magnitude as the temperature rises. In the following section we 
shall give a thermodynamic argument to support the idea that e and 5 should in fact be chosen 
to ensure that <§ < as T -> 0. We shall see that this can indeed be done, provided that they 
are greater than some critical value. 

Before proceeding, we add one further comment. The attentive reader might object to our 
previous remark about the mass being analogous to an order parameter, on the grounds that 
the Coleman-Mermin- Wagner theorem p3| , p4| applies at finite temperature for a system in two 
spatial dimensions. This theorem forbids the existence of an order parameter which breaks a 
continuous symmetry of the system. However, the Kosterlitz-Thouless mechanism allows the 
global T3-symmetry of QED 3 to remain unbroken despite the dynamical generation of mass. 
According to Witten ||25| , the existence of dynamically generated fermion mass, which is the 
usual physical consequence of a non-zero order parameter, shows that even though (iptp) is not 
quite an order parameter, it is very close to being one. 
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4 Numerical Results 



In this section we shall solve numerically equation (|TTJ) at N = 1 for various values of 5 and e. 
We shall see that for the specific values S = S c and e = e c , ^ becomes negative for small T/a. 
As we have stated in the previous section, we believe our results to be unphysical below 8 C and 
e c ; in the next section we shall support this view by a simple thermodynamic argument. The 
behaviour of E as a function of T is insensitive to small variations of N |T7], |TJ[] . The larger the 
N, the smaller the T r and the more difficult it is to see the numerical solutions. 



In our numerical evaluation of (11), we must consider the T = case as well as that of T ^ 0. 
To solve for s(T — 0) we use an equation which was derived in |Tj| for the zero temperature 
case of @, retaining only the longitudinal part: 



Q(oo, s) __ 2 
2nN = (2it) 2 N 



(0.125) 2 \ 7i 0.125 

In 



(23) 



(0.125) 2 + s 2 J 2s (0.125) 2 + s 2 V0.125, 

Before using this equation, we need to replace iV by N/2. This is because ( p3|) was only derived 
for the longitudinal part of @. Since at T = the transverse and the longitudinal contributions 
are identical we need only double the longitudinal contribution. 

Evaluating (^) numerically yields a value of s(T = 0) = s ^~°- > = 0.0578; note that this 
result as well as (|23|) itself are cutoff independent, due to the fact that both cutoffs go to 
zero at T = 0, where there is no infrared divergence. There is one more crucial property of 
(|2"3|); lim T ^ Qo( a ) = oo, where Qo(a) = Q(a,s = 0). Combining this result with some other 
considerations in Appendix B, will lead to conclusions about the critical temperature and about 
some properties of a solution s = s(T) of the S-D equation. 

To solve (|Tl"D numerically at T ^ for both 5 and e dependence, we choose a value of S or 
e and seek to determine l/a c = To do this we know from Appendix B that Qo(cl c ) = 27r 

at N = 1. Therefore we start by varying Qq as a function of a until the condition Qq = 2tt is 
satisfied, the value of a at that point being a c . Once we have determined a c , we start looking 
at values of a > a c for which Q > 2ix . For Q < 2n we know from Appendix B that E = 0. 
At each value of a we consider, we vary Q(a, s) as a function of s until ([11]) is satisfied, namely 
such that Q(a, s) = 2tt for N=l. Using each value of s = E/a we evaluate for a given a, as well 
as for the T = solution, we can construct full temperature dependent solutions of flTT|) for a 
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particular value of e or 5. 

The first set of numerical solutions are those for the integral cutoff method, as shown in 
figfj. For each of these solutions we give a value of and we calculate the ratio r = 2 ^^~°' > ■ 
As one can see for all values of e, except for e = 2.5 and e = 2, the gradient ^ is positive 
for small T/a. As e decreases we see that this effect becomes more pronounced. We also 
see that as e decreases, r decreases (as shown in table [I]); this is due to the increase in the 
critical temperature, while at the same time still retaining a cutoff independent T — > limit. 
Both these effect are consistent with the solutions of the previous section for extremely small 
e, where we have found that T c ~ with h = In (~\ Also we see for these extremely small 
e solutions, that ^ — > oo as T — > (see fig.llla). 

The second set of numerical solutions are those using a mass cutoff and are shown in fig.[| 
We give values of r and ksT c for each solution of a specific 5 in table [| It is important to 
notice that again there is a transition: below a certain value of 5 we see that ^ is positive for 
small T/a. A positive gradient case is shown for 5 = 0.05. Again we notice that the r- values 
fall with decreasing 5, consistent with the extremely small S solution of the previous section. 
One should notice that for the mass cutoff, if we look only at r values for solutions which do 
not have a positive gradient, we see that r changes slowly, namely by a factor of order 0.7 over 
an order of magnitude in 5 (5 ranging from 2 to 0.25). We see that r = 6.2 for the value of 
5 = 0.5, which is very close to the value of r given in |n| of 6.17. This is because S is at a 
value close to the corresponding longitudinal plasmon mass of ^-^a, so making the transverse 



contribution roughly the same as the longitudinal, as in [IS 



5 A thermodynamic argument for ruling out 5 < S c and 
e < e c solutions 

We now look for a simple thermodynamic argument to rule out solutions with positive gradient 
at small T/a. We consider a gas of fermions, which we take to be dilute, so that we may ignore 
any interaction term. The chemical potential of a dilute gas is negligible and may be discarded. 
We shall also assume the mass to be temperature dependent. The partition function takes the 
usual form for a non- interacting gas of fermions (with fj, — 0), except that now cu = (p 2 + E 2 ) 1 / 2 
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C 


9 5 


9 


1 ^ 


i 


^ 


9 s ) 


T 

J- c 


016 




023 


026 


035 


0445 


r 


7.23 


6.42 


5.03 


4.45 


3.30 


2.60 



Table 1: The values of r = 2 ^ T „ ^ for different values of the integral cutoff e. E is the 

2 

dynamically generated mass and T c the critical temperature measured in units f^. For the 
allowed region of values e > e c , e c ~ 2, the average value of r is about 6.8 in agreement with 
previous papers. 



5 


2 


1 


0.5 


0.25 


0.05 


T c 


0.015 


0.0164 


0.0185 


0.0215 


0.031 


r 


7.70 


7.05 


6.25 


5.38 


3.73 



Table 2: The values of r = 2S fc ( J~ 0) for different values of the mass cutoff 5. £ is the dynamically 

2 

generated mass and T c the critical temperature measured in units j^. For the allowed region of 
values 5 > 5 C , 5 C ~ 0.2, the average value of r is about 6.6 in agreement with previous papers. 
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is temperature dependent: 



\nZ= 2VN J J4jL \j3u(T) + 2 In (l + e"^ (T) )] (24) 

where 1/ is the volume of the system and N is the number of flavours. We shall calculate the 
entropy of this gas. It is convenient to divide In Z into two contributions, In Z = In Zq + In Ze, 
where In Zq is the contribution from the ground state and In Ze is the contribution from the 
excited states. We first consider the ground state contribution. We calculate the frequency 
integral imposing an ultraviolet frequency cutoff, Ac, since the integral in InZo is ultraviolet 
divergent. We then have 



1 7 VN ? 

In Z G = 



7T 

,3 



A 3 G £ 3 (T) 



(25) 



Aq should be thought as an (infinite) constant arising from the zero-point energy contribution 
at T = 0; since it is constant it will not affect the entropy. We know that the entropy is related 
to the partition function by S = 9T q^ z (in these calculations we shall set = 1). From this 
we find that the entropy of the ground state takes the form: 

VN „, N dE 

Sg = £ 2 TV- 26 

7T dl 

We now look at the contribution arising from excited fermion states. To perform the integral 
analytically we expand out the logarithm in (|24|) for /3E(T) > 1 retaining only the first term. 
This is allowed, because we are interested only for T/a small and E(T = 0)/« ^ 0. Thus we 
arrive at the following form for the In Ze term: 



AVT 2 N /£ 
2tt \f 

It follows that the corresponding entropy, Se , at small T/a is given by 



lnZE = 2L i^{f + 1 ^ e ~^- (27) 



_ 4VNT 2 /3S(T) S 2 (T) E dX\ ^ 

For very large E/T we notice that (|^) is suppressed by a factor of exp(— E/T), so at large 
E(T)/T this contributes little to the entropy of the system. It is Sg which is the dominant 
factor. We therefore conclude that, although our model may be an over-simplification of a 
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real physical system in QED3, positive ^jp at small T/a leads to negative entropy, which is 
clearly unphysical. The consequence of this is that there is a boundary between physical and 
non-physical solutions in the parameter spaces of the two types of solution. From this we may 
conclude that the contribution from higher orders in 1/N, if indeed it regulates the infrared 
divergence, must be above a certain magnitude. 

There is one point that needs to be addressed before we finish this section. If ^p 7^ at 
T = 0, fl26f ) implies that S 7^ at T = 0, which is unphysical: when T — ► the 3rd law states 
that S — ■> 0. From our numerical data we strongly suspect that for all values of S and e (except 
5 = 0, e = 0) ^p = at T=0, when all modes are considered. Although we have been unable to 
prove this, we are able to present an argument showing that ^p ^ is independent of e and 5 in 
Appendix B, so that ^p is the same in all cases (except when e(6) = ). We show in fig.f| 
what we think happens as T — > 0, where the gradient of all the curves goes to zero as T — > in 
accordance with the 3rd law. Although this behaviour is likely to be true for very small T/a, 
^jp is significantly positive at small enough T/a for e < e c and 5 < 5 C , so that e c and S c are 
good estimates of the boundary between physical and non-physical regions; the contribution 
from fl2§D is very small, when compared to (|26f) in the region where we see positive ^p. 



6 An approximate form for E and estimates of 5 C and e c 

In this section we obtain an analytic representation of S as a function of T and e(5). This 
will give us a simple way of estimating e c (5 C ) and it may also be useful in applications to 
phenomenology, for example to a calculation of the temperature dependence of the Meissner 
effect (although we shall see later that a mass cutoff is not consistent with the U(l) x U(l) 
model of superconductivity, []14[ ). From our numerical solutions (fig.0 and fig.|3] ) we see that 
S as a function of T looks very much like the section of an ellipse, centered about some value 
T , in the quadrant S > and T > 0. Such an ellipse takes the form (see fig.§|a, fig.03) 

m 2 (T) = 1 + 2T ° T — (29) 

' T[ 2 -2T T^ T' C 2 -2T T' C 1 ; 

in terms of the two parameters T and T' c and where m(T) = s pp jL ■ Fitting (p9|) to the 

numerical data, we have found that it is better to use both T an d T' c to fit S(T), instead of 
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fixing the value of T' c at the calculated T c for a particular e (5) and varying only To. To indicate 
the accuracy of the fit we show in fig.|]a and fig.|]b both the approximation and the numerical 
data for one chosen value of e and one of 5. Fitting both T and T' c we are able to build up 
graphs of T for both cutoffs (see fig.[7]a and fig.[7|b). Using the simple condition that T < 
implies 4iS < 0, we are able to estimate e c and <5 C ; we find e c ~ 2 and 5 C ~ 0.2. 



7 Discussion 

Although in QED3 both types of cutoff lead to roughly the same behaviour for the mass gap, 
as a function of T, the mass cutoff is ruled out, if applied to the U(l) x £7(1) model of high- 



ly superconductivity [Tj]. The problem is that the transverse part of A Mi , (the statistical 
gauge field propagator) in this model must be massless. If A M „ has a transverse mass, it 
no longer generates (via an analogue of the Higgs mechanism) a transverse mass in D^, the 
electromagnetic propagator and there is therefore no Meissner effect. This also suggests that 
infrared regularizing contributions from beyond leading order in 1/N may come from the vertex 
function in (§), if QED3 is indeed a model of high-T c superconductivity. 

Although we found (Appendix A) that the transverse contribution to the photon propagator 
still behaves at pob = for small |p| as 1/p 2 , even if massive fermions are used in the vacuum 
polarization, our expressions for IT^, remain useful. It will be interesting to see what effect 
a photon propagator involving 11^, (instead of 11^) has on the calculation of a 3-momentum 
independent mass, S. By looking at (A. 21), we see that as /3E — > 00, the plasmon mass tends 
to zero in the longitudinal part of the photon propagator at p fc = 0. It may be that instead 
of monotonically decreasing with /3E, at large enough /3E, the contribution from the zeroth 
longitudinal mode to the gap equation is enhanced. If this behaviour is present in the zeroth 



longitudinal mode, it might cause (|TT| ) to favour smaller values of S than would be expected, 
if, as in the case of IL^, the zeroth mode plasmon mass contribution was independent of E. 
It might also be interesting in the context of the infrared problem to look at the effects of 



the wavefunction renormalization given in |T7|, where Ai, the physical mass, is related to E by 



»™ - life < 30) 

We note that the integrand in the equation for A(j>f) has a pole at qb = pj — kf = 0, where 

17 



kf is the loop momentum (see equation (6) of |T7|]). We strongly suspect that this will lead 
effectively to an infrared divergence, which will then dominate the equation for A(pf) for very 
small cutoff. For the case of a 3-momentum independent Ai, one should be able to derive 
an expression relating A(pf) to M.. From this one could see what effect A(pf) has on A4 for 
small cutoffs, the hope being that M. may be finite for T / without a cutoff. One should 
note, however, that we may have to enforce the Ward-identities at finite temperature, which 



the wavefunction renormalization approach in [17]] does not do, in order to regulate the infrared 
divergence. 
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Appendix A - Calculation of the polarization tensor with a fermion mass to 

leading order in 1/N. 

In [[OJ it was found, for the polarization tensor IT^ calculated for massless fermions, that 
there was an infrared divergence in the Schwinger-Dyson equation for the fermion mass. This 
infrared divergence comes from the zeroth frequency part of the photon propagator which goes 
like 1/p 2 for small momenta at T ^ 0. In this appendix we shall calculate the polarization 
tensor IT^, with fermion propagators of mass S; we use the superscript £ to distinguish this 
polarization tensor from that in the massless case. Again this calculation is to be done at 
finite T and in the Matsubara formalism. It is precisely the low-|p| behaviour of the resultant 
photon propagator, A^ v at pob = 0, that we shall be interested in; we want to see if the infrared 
divergence is regulated by the replacement of the massless fermion propagators with massive 
fermion propagators in our calculation of iT^,. 

To calculate IT^, to leading order in 1/N we need only consider the contribution from the 
diagram shown in fig.Al. 
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From fig.Al we deduce that the polarization tensor takes the following form 



i-rv i \ a f d 2 k ( 1 1 \ 

n -W = fl?y^ ta ^'S^+Wi^"SFwJ- (A ' 1) 

In our calculation of II?, we take £ to be constant with respect to 3-momentum. If E 
was to vary with respect to frequency or momentum, we would at this point need to insert 
a variational ansatz and solve integral equations for the variational parameters derived from 
the Schwinger-Dyson equation for the mass, S; or we would have to treat (A.l) as a separate 
integral equation. By use of Feynman parameterization and other manipulations we arrive at 
the following expression for IF^,: 

s _ 4a r 1 ^ r d 2 l 2y„ 6^ 



Jo 



Uflu (3 Jo ^ J (2vr) 2 (P + p 2 b x(l -x) + £ 2 ) 2 (/ 2 + p 2 x(l - x) + £ 2 ) 
2(x(l - x))(p 2 b 5^ - PfAPvh) (1 - 2x)(l fl p 1/b + l u p^ b - (pb.^S^) 
{l 2 + p 2 b x{l-x) + T 2 ) 2 {l 2 + p 2 x{l-x) + T 2 ) 2 1 ' j 

where / is the shifted 3-momentum, related to kf by I = kf + xp b . For our purposes we need 

only consider II? and IIq . 

Before proceeding any further there are some immediate points that need to be raised. The 
first is that, if we were to take the limit T — > of (A.2), the first two terms would cancel after 
regularization. This is because the first two terms have the same tensor structure, namely both 
are proportional to 5^, and using dimensional regularization one is able to show that these two 
terms are equal in magnitude and opposite in sign. But now consider the case when T ^ 0: 
in this case we have a sum over k f components in l instead of an / -integral and so there is 
a preferred direction in 3-momentum space. The upshot of this is that the dimension of the 
/-integrations is effectively reduced and the first term is no longer proportional to 5^ v . 

In the case of IT? the first two terms in (A.2) do cancel at T ^ 0, due to the /-integrals 
having the same form as those for T = 0, although now two dimensional instead of three 
dimensional. For the IIq terms the situation is slightly trickier, for the first term of (A.2) 
instead of having Ijlj as a factor has /q. It is this term in ITq which causes the first term to 
be no longer proportional to 5^, and this means that for IIq the first two terms in (A.2) no 
longer cancel. 

At this point we should ask ourselves the question: is there any constant, leading order term 
when we set p b = 0? If we look at II? , the answer is no, for when we set p b = only the first 
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two terms survive, which cancel. For I1q , the answer is yes, for although only the first two 
terms survive, they do not cancel. This is precisely why, as we shall see in our calculation of 
11^, there is a plasmon mass only in the longitudinal part of the photon propagator at p ob = 0, 
p = 0. Exactly the same argument holds for U^, the massless case. Although this is true, we 
are still interested to see what is the actual low-|p| behaviour the p ob — part of the photon 
propagator. 

Continuing our analysis, we write both ITq and as: 

!F -— Cdx ( — 
11 00 - p J Q dxJ (2?r)2l 

and 

d 2 l 



■\Sl + 2(1 2 + £ 2 + x(l - x)p 2 b )S^ + (1 - 2x) Pob S* 2 



(A.3) 



where 



J Q dxJ j^[2x(l-x)( P 2 b 5 ij - PiPj )Sf + (l-2x) Pob Sf - 5^(1 -2x) Pob Sf] (A.4) 



s?= E 

n=—c 

oo 

sf= E 



i 



(z 2 + i 2 + s 2 + x(\-x) P 2 y 

^0 

(z 2 + i 2 + s 2 + x(i-x) P 2 y 



(A.5) 



From (A.5) one can easily deduce the following relationships 

1 dSf 1 



1 dSf 
2y dy 

s _ _pds? 

b2 ~ 2 doo 



dSf 



2|1| 



2x(l - x)p 6 dp b 



(A.6) 
(A.7) 



where y 2 = l 2 + X 2 + x(l — x)p 2 and a; = n(2mx + 1). Using (A.6) and (A.7) we are able to 



express 11™ and Uf • 



d\l\ 



n 



(2*) 
/3 Jo 7 (2tt) 



|1|S' 1 S + (|i| 2 + s 



2 ,0Sf , p 2 dS\ 



+ 



|1| - (l-2x) 



pp ob dS 1 



d\\\ p dp ' 2 duj 



-(P6^--W)g^|l| + (i _ 2x)^^5ij 
p dp 2 duo 



We can evaluate the sum Sf by the usual complex integration, so obtaining the result 

s /3. u + ipy uj-i/3y 
S 1 =4fil<x*(—j— ) + cot(^— )] 



(A.8) 
(A.9) 

(A.10) 



21 



where 



Using (A. 8), (A. 9) and (A. 10), as well as regularizing the integrals so that the ultraviolet 
singularities cancel, one is able to evaluate the |l|-integral. This gives us final expressions for 
Hoo and ng, 

n s = nf - 4nf - nf + n 4 s (A.n) 
Pb 

n§ = nf (dy - ?f ) + n^, (A.12) 
Pb 

1 x(l - x) sinh(/j(p^(l -x) + E 2 ) 1 / 2 ) 

X (a;(l-a;)p 2 + S 2 )V2 D^nJ) 
Z" 1 / n , sin(2xm7r) 

/ l-2xrfx-^ ^ A.13 

— / dxln(4D*(x, Pb ,f3)) 

TTfJ JO 

1 dxsmh(f3(p 2 b x(l -x) + S 2 ) 1 / 2 ) 
(pl(l-x) + ^y/Wl(x,p bj (3) 

Dl{x lPbl P) = cosh 2 ^(pMI-^ + S 2 ) 1 / 2 ^ _ sin 2 nmx (A W) 

We now examine whether the expressions (A. 11), (A.12), (A.13) and (A.14) give the correct 
limits. We can easily show that lim s _+ong = and lim s ^ riQ = n 00 , where and n 00 , 
the results for the massless calculation, are given in |14, When E — ► oo we require that 
II^ — > and IIoo — > 0; by looking at (A.l), we can see that this requirement must be true. 



n? = 


apt r 
2tt Jo 


n 2 s = 


( am\ 
I 2/3 J 


n 3 — 


a r 1 
7T/3 Jo 


114 — 


a/3E 


2tt 



and 



Also to be consistent with the results at zero temperature these terms must vanish in the 
E — > 00 limit (although this does not rule out the presence of a term independent of E and 
proportional to a/ (3 which vanishes in this limit). For /3E 3> /3pb and /?£ > 1 we have the 
following: 



ni^^- n 2 ^e-'^pntO (A.15) 

12L7T p7r 

*-?(!) n < -f (!)■ 

From (A.15) it is easy to show that ITq and IK- have the correct E — > 00 limit. 

We now need to deduce the low |p|-behaviour for the photon propagator. In the Landau 
gauge we know that the photon propagator must take the form: 

v pl + Til{p h ) + rf + ng(p 6 ) lA ' iDJ 
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where A^ u and B^, the longitudinal and transverse projection operators, are given, to leading 
order in 1/N, by 



V=^-^)(g)(*.-^j (A.17) 
B^ = S^U ii - 1 ^j8 jv . (A. 18) 

11^ and 11^ are related to Hfj and Uq Q by the following formulae (for further discussion 



sec 



14,1 



n5 = n0 nl = n|-n^| (A.19) 

Using (A. 19) we can deduce the zero frequency forms of 11^ and II |j in terms of Ilf , Ilg and 
Ilf . Noting that at p b, nf = for all |p|, 

n^ ob = o) = nf + nf u%( Pob = o) = uf (A.20) 

From (A. 11), (A. 12), (A. 13) and (A.20) we find that to leading order in |p|: 

*i = — n rn, ^ a tt^t + I^^tt^ (A.21) 



Note that unlike (A. 15), we have not taken /3E 3> 1, but rather /3S S> /3|p| at po& — 0; this 
gives us factors of tanh/3S and ln(2 cosh/?S/2). From (A.21) one can plainly see that A^ u 
goes as 1/p 2 in the transverse part for small |p|. This is exactly the same behaviour as A^, 
calculated with massless fermion propagators; therefore the inclusion of the finite fermion mass 
in the vacuum polarization does not remove the infrared divergence in the gap equation. 

Appendix B - Analytic properties of the solution to the gap equation. 

In this appendix we are interested in those properties of the solution S = S(T) of the 



S-D equation in the form (|TT| ) which can be found analytically. In the analysis that follows, 
calculations will be given for the integral cutoff case, but the properties found will be true of 
the mass cutoff case also. For any finite number of flavours N, we shall prove the following 
points: 

(a) There is a non-zero critical temperature T c . 
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(b) T c is unique. 

(c) For T > T c , s is identically zero. 

(d) At T = T c , § -> -oo. 

(e) At T = 0, is cutoff independent. 

It should be mentioned in passing that all these proofs apply also to the solution for constant 
S, where only the longitudinal part of ( |TT| ) is retained as in ||19| . 

We define Qo(T) = Q(T, s 2 = 0). For the purposes of this appendix we first need to prove 
that Hag* < and ^1 < . 

To show -j^r < 0, instead of using equation (pTT| ) it is much easier to go back to (y), ( JT0| ) 
and define Qo(T) as: 



oo „ 

Qo(T) = aJ2] 

m=l e 



2xdx 



X 



x 2 + (27rm) 2 + 0.125a[x 2 + (27rm) 2 ] 1 / 2 
1 1 



x 2 + [2vr(m+ 1/2)] 2 x 2 + [27r(m - 1/2)] 5 



+ 



00 ax<ix 



a; 2 + tc 2 

Differentiating Q with respect to a yields 



+ 



1 



x 2 + (3 2 U 1 x 2 + (3 2 U 3 



(B.l) 



dQo 
da 



E 

m=l 



2x[x 2 + (27rm) 2 ]dx 



x 



(x 2 + (2vrm) 2 + 0.125a[x 2 + (2™) 2 ] 1/2 ) 2 
1 1 



+ 



x 2 + [2n(m + 1/2)] 2 x 2 + [27r(m - 1/2)] 



+ 



00 x 3 <ix 

X 2 + 7T 2 



1 



+ 



1 



[x 2 + /? 2 n^] 2 [x 2 + p 2 u^ 2 



01 



(B.2) 



It is now obvious that, since x > e > 0, -^p > and consequently = — ^ < for any 
T / 0. Noticing that the quantities /3 2 ITq, /3 2 I1q are always positive for x > and a > 0, we 
see that 



dQ(T,s 2 
ds 2 



x 



-« 3 E 

lm=l 
1 



2xdx 



x 2 + (2vrm) 2 + 0.125a[x 2 + (2vrm) 2 ] 1 /2 

1 



+ 



[x 2 + [27r(m + i)] 2 + a 2 s 2 ] 2 [x 2 + [27r(m - ±)] 2 + a 2 s 2 ] 2 
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+ 



00 xdx 



1 1 



x 2 + (3 2 Ii l + x 2 + (3 2 Ul 



(B-3) 



[x 2 + 7T 2 + a 2 s 2 } 2 

is always negative. There are two more properties of Qo(T) which will come in use in the 
following, namely lim^o Qo (T) = oo which we discussed in section [| and lim^oo Qo{T) = 
which can be easily seen from (|B.1| ). We are now ready to proceed to the discusssion of points 
(a) to (e). 

(a) Consider the surface Q(T, s 2 ). A critical temperature T c is a value of T such that s = 
and at the same time the S-D equation is satisfied, that is 2ttN = Q (T). In geometrical terms 
it is the T-coordinate of a point which belongs both to the surface Q(T,s 2 ) and to the line 
which is defined by the intersection of the planes s 2 = and Q = 2nN . We have proved that 
the curve Qo(T) is decreasing with T and lim^o Qo(T) = oo, lim^oo Qo(T) = 0. Since Qo(T) 
is a continuous function of T, it follows that for any given N there is a critical temperature 
T c = T C (N). For any finite N, T c ^ 0. The larger the iV is, the smaller the T c is in agreement 
with ||,@. 

(b) For any T > 0, Qo(T) is monotonically decreasing with T. Thus the S-D equation 
2tiN = Qo(T) defines a unique T c . 

(c) For any finite non-zero T and any real non-zero s, Q(T, s 2 ) < Qo(T), because Q(T, s 2 ) is 
a monotonically decreasing function of s 2 . For T > T c , Qq{T) < Qo(T c ) = 2ttN, because Qo(T) 
is a monotonically decreasing function of T. Therefore Q(T, s 2 ) < 2tcN which means that for 
T > T c there is no real non-zero s which satisfies the S-D equations. Combining (c) with (b), 
we see that, for T > T c , only complex solutions for s can exist. This would lead to unstable 
particles or tachyons. However above T c we can choose another solution to ([5|), namely s = 
identically. Thus QED 3 reproduces the superconducting feature of the gap closing above some 
critical temperature. 

(d) Consider that the curve s = s(T) is a solution of the S-D equation, that is 2ttN = 
Q(T, s(T) 2 ). It is s(T c ) =0. In a neighboroughhood of T c we can make the approximation 



Q(T lS (TY) = Q (T c ) + 



ds 2 

Rt + Rs 1T 



(T-T c ), (B.4) 



where R T = dQ{ If 2) , < and R s = dQ ^f ] , < 0. Recalling that Q (T t 



T=T c ,s 2 =0 



T=T c ,s 2 =0 
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2irN } we can rewrite ( |B.4p as 

2s(T c ) ^ 
v ; dT 



(B.5) 



The right hand side of the above equation is a negative finite number for any finite temperature. 



= — oo. 



Thus the fact that s(T c ) = + implies that 4l 

(e) We now turn our attention to arguing that ^ is cutoff independent. For this purpose 
we shall approximate ITq and I1q in two regions. In the first region we approximate both I1q 
and I1q by their low-|k| behaviour, namely /3 2 ITq ~ a fx 2 and (3 2 Uq ~ 2 ln2a/7r; this region is to 
be defined by k-momenta whose magnitude is below A. In the second region we approximate 
both /3 2 ]Tq and /3 2 ITq by ax/8, their behaviour at high-|k|; we define this region by |k| > A. To 
make both ITq and IIq continuous, A should be chosen to be 16ln2/n. In this argument we shall 
be considering e < A, although it should be straightforward to extend our analysis to e > A 
(care being taken to include e in other modes for e > 2n). Using these approximations we are 
able to write the contribution from the full zeroth mode as 

' [F (e, A, a, s) + G (A, a, s)] (B.6) 



2nN 
where 



, . ax dx I 1 1 \ , . 

o(£ ' ' a ' S) = X x 2 + a 2 S 2 + 7r 2 \x 2 + afx 2 + x 2 + 2aln2/J ( } 



and 



Go(A ' a ' S) = 2 / A x 2 + 0125ax C^ + aV + vr 2 '- (B " 8) 
The integrals in these expressions can be evaluated analytically 
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a/7t 2 + a 2 s 2 



(B.9) 
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Since we are interested in the low T behaviour of both Fq and Go, we expand up to second 
order in (1/a) for both terms, replacing leading order terms in Go by 1(0, 0.125a, as). 



7T 



4s 2 In 2 



F (e,A,a,s) ~ ^F{e,A,s) = — 
a z a 1 

G (A, a, s) ~ -^G(A, s) + 2a/(0, 0.125a, as) = 2aJ(0, 0.125a, as) 
a 

A 0.125AN 



a 2 (s 2 + (0.125) 2 ) VO.125 s 
We define F c as follows 



(B.10) 



fB.ll^ 



]T 2al(\2irm\, 0.125a, (aV + [2vr(m + 1/2)) 2 ] 1/2 ) = Q(oo, s(T)) - F c (a, s) (B.12) 



m=— oo 



where Q(oo, s(T)) coincides with the formula of section [|, if s(0) is replaced by s(T), so that 
Fc(oo, s) = 0. Now we can write (|TT|) as 



2nN = ^[F(e, A, s) + G(A, s)) + Q(oo, s) - F c (a, s). 



(B.13) 



Differentiating the above equation with respect to T yields an equation for at low T/a 



= — [F(e,A,s) + G(A,s)] + — + 



a" 



ds 
dT 



-F c (a,s 



ds ds 

dQ(oo, s(T)) ds 



ds 



dT 



(B.14) 



where F' c = Since ^j, ^ are finite at s > 0, e > 0, this means that all the cutoff 

dependent terms vanish, for s is non-singular at T = and we get a new gap equation which 



is cutoff independent. This means that is independent of e. To show that J| 

one must show that F' c (oo, s) = . Then one is left with the gap equation 



T=0 











ds dQ(oo,s) 
dT 



ds 



(B.15) 



for which the solution is = for finite s. 
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Figure 1: (b) The extreme solution for the mass cutoff case (h = 500). 
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Figure 2: Numerical solutions for the integral cutoff e. The curves show the scaled dynamically 
generated fermion mass as a function of the scaled temperature for various values of e and for 
N=l. 
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Figure 3: Numerical solutions for the mass cutoff 5. The curves show the scaled dynamically 
generated fermion mass as a function of the scaled temperature for various values of S and for 
N =1. 
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Figure 4: Each curve corresponds to a different value of e (or 5). All the curves have zero 
gradient as T — > 0, liniT_^o 3^ = 0, satisfying thus the third law of thermodynamics. 
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Figure 6: (a): The diagram shows the fitting of the ellipse approximation (j^J = 0.00334233 
0.00546513 (J) - 12.8907 (J) 2 to the points we have found for e = 2.5. 
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I) -0.00334233 



Figure 6: (b): The diagram shows the fitting of the ellipse approximation f^J = 0.00334233 — 
0.0497448 (J) - 7.05771 {^f to the points we have found for 5 = 0.5. 
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